Supermassive Black Holes in Galactic Nuclei 
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ABSTRACT 

We discuss the link between the observations of distant quasars and those of massive dark 
objects in the cores of many local galaxies. We show how the formation of early black holes 
gives rise to the luminosity function of high z quasars, while it imprints into their dark local 
relics a related shape of the mass-dispersion correlation. We propose that in its lower section 
the correlation slope will tell the (otherwise uncertain) strength of the feedback effects from the 
quasar radiation on the host galaxies. 
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Introduction 



Black hole physics - galaxies: active - galaxies: interactions - galaxies: nuclei 



Two news have recently kindled the field of the 
quasars and active galactic nuclei. 

The first concerns the farthest objects. Not 
only single quasars (QSs) have been detected at 
redshifts out to z = 6.28, but also the statistics at 
z s» 5 has improved to the point of outlining the 
bright section of the luminosity function (Fan et 
al. 2001a, 2001b). 

These findings confirm that the population of 
the optically selected objects goes through the 
most sharp and non-monotonic of evolutions. The 
comoving density of the bright sources rises on 
the scale of a few Gyrs from the Bang, to peak 
at around z » 3 (Shaver et al. 1996); later on 
it turns over and falls by factors 10 -2 toward us 
(Boyle et al. 2000). A similar message comes from 
the radio band, see Jackson & Wall (1999). 

The second piece of news concerns the cores of 
many local galaxies, where massive dark objects 
(MDOs) ranging from a few 10 6 to a some 10 9 Mq 
had been detected within regions from a few to 
tens of pes, see Richstone et al. (1998). Now (Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000) such 
masses are found to correlate tightly with the ve- 



locity dispersion in the body of the surrounding 
host galaxies. 

The MDO masses fit into the framework 
provided by the long standing arguments (see 
Lynden-Bell 1969) that indicate gravitational con- 
traction rather than thermonuclear burning to be 
the dominant source of QS output. We just recall 
this conclusion to hinge upon the high bolometric 
power L > 10 45 erg s _1 of many such sources, 
and on the high compactness with sizes down to 
R ~ 10 15 cm indicated for some of them; it also 
requires an overall efficiency up to rj ~ 10 _1 for 
conversion of gravitational into radiative power. 

But then the argument may be carried on to 
evaluating the masses involved (Cavaliere et al. 
1983); in terms of M s = M/1O 8 M these read 



M fi 



3(L 45 At-!Rvi/v-i) 



1/2 



(1) 



where we have used L45 = L/10 45 erg s _1 , 
i?i5 = i?/10 15 cm, and the source life time 
At = At-i 10 _1 Gyr. Such masses are consis- 
tent with the MDO observations. 

Here we intend to link the largest values of M 
with the early QSs, and to discuss how such a 
relation is to be extended into the range of lower 
redshifts and smaller masses. 
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2. Black holes and their environment 

The accreting black hole (BH) paradigm not 
only best accounts for the small sizes and the top 
efficiencies required, but also provides to the spent 
masses the stability of a terminal configuration 
(Rees 1984). Nailing down the indication from eq. 
(1), a BH keeps full record of the mass 

M = J dtL^/rjc 2 (2) 

accreted over the life of an active galactic nucleus. 
Meanwhile, the induced bolometric luminosities 

LKitj^Am/At (3) 

are tuned (even at constant rj <~ 10 _1 ) by the mass 
Am accreted over a time At, and cover a wide 
dynamic range: from Eddington self-limiting con- 
ditions governed by radiation pressure that yield 
L ~ Le ~ 10 46 M 8 erg s^ 1 , to supply-limited ac- 
cretion (Cavaliere & Padovani 1989) that easily 
allows only sub-Eddington emission. 

Thus strong gravity is not enough. Equally 
important are the environmental conditions in, 
or surrounding the host galaxy; these can drive 
widely ranging accretion rates from 1CP 3 or less, 
up to some 10 2 M Q yr _1 . In addition, it is the 
cosmological change of the environment that qual- 
ifies to govern the QS evolution; in fact, both oc- 
cur on time scales t ev <~ a few Gyrs such that 
rjtE -C t ev < Hq 1 holds. 

In turn, the environmental conditions are de- 
scribed by the other paradigm, the hierarchical 
growth and clustering of dark matter (DM) ha- 
los, wherein the galaxies constitute lighter bary- 
onic cores (White & Rees 1978). This implies sub- 
stantial dynamical events to occur, at early z in 
the strong form of merging between comparable 
subgalactic units, later on as milder interactions 
of galaxies in groups. 

All these dynamical events tend to break on 
scales of kpcs the axial symmetry of the galac- 
tic gravitational potential, or enhance its steady 
asymmetry; relatedly, the specific angular momen- 
tum j providing support to the gas in the central 
kpc of the host is not conserved, rather it is trans- 
ferred to the massive DM component. Thus the 
necessary condition is provided for destabilizing 
and funneling inward a sizeable gas fraction. At 



smaller scales dissipative processes take over to re- 
distribute j (Haehnelt & Rees 1993), and cause the 
gas to reach the nuclear accretion disk and grow 
new BHs or refuel the old ones. 

To link BH growth and QS luminosities we fol- 
low Cavaliere & Vittorini (2000) and disentan- 
gle the triggering dynamical events into two main 
regimes, roughly divided by the epoch of group 
formation zq ~ 2.5 ± 0.5, depending on cosmolog- 
ical/cosmogonical parameters. 

3. QSs in forming spheroids for z > 3 

The self-limited regime occurs mainly at epochs 
before zq, when galactic spheroids are built up 
through major merging events between halos with 
M/j < 10 13 M Q . These events destabilize large 
amounts of gas, but also replenish the host struc- 
tures with fresh supplies and sustain the gas 
amount at the cosmic level m « 10" 1 M^. 

As a consequence, central BHs can form and/or 
accrete rapidly, growing by Am ~ M over dynam- 
ical times close to the Salpeter scale, td yn ~ f]kE\ 
so after eq. (3) they attain Eddington luminosities 
L ~ L E cx M. 

In turn, M is related to M^; two specific mod- 
els bracket the processes involved, see Cavaliere & 
Vittorini 1998 (CV98), Hosokawa et al. (2001). 
Haehnelt & Rees (1993) considered BH coales- 
cence in parallel with merging of their halos; this 
process (hereafter model A) is described by the 
simple scaling 

M « 10~ 4 M h . (4A) 

Alternatively, Haehnelt, Natarajan & Rees 
(1998) proposed the feedback-constrained model 
B where the scaling reads 

M^{l + zf' 2 Ml%. (4S) 

This is because during halo merging a central BH 
may also accrete gas up to the limit 

e Lstdyn ~ G M h m/r ; (5) 

this is set by gas unbinding from the halo poten- 
tial well, due to the deposition of an (uncertain) 
fraction e ~ 10~ 2 of the QS radiation (Silk & Rees 
1998). 

With both models (4A) and (4B) the early QSs 
are expected to grow in average luminosity and in 
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number, tracking the development of protogalac- 
tic halos over the range from Mh ~ 10 10 toward 
1O 13 M . The evolving halo mass distribution is 
widely taken in the form Nps(Mh, z) first pro- 
posed by Press & Schechtcr (1974); the positive 
term of its time derivative provides the rate of 
halo formation, and yields (see CV98) the lumi- 
nosity function (LF) in the form 

N{L, z) dL = At d t N PS (M hl z) dM h , (6) 

with the prefactor At w rjtE accounting for the 
limited source lifetime. 

Fig. 1 shows the optical LFs provided by such 
models. We stress that model (4B) by the very 
means of its non-linear transformation stretches 
the halo distribution into flatter, more fitting LFs 
and predicts a stronger (negative) evolution. By 
the same token, it also associates bright QSs with 
the largest galactic halos, consistent with the data 
discussed by Hamilton, Casertano & Turnshck 
(2000). 

Both models are normalized to the data at 
z ~ 4. But model (4B) privileges the upper halo 
range where N(M h ) decreases steeply, so it pro- 
vides BH numbers smaller and naturally close to 
one large BH per actively star forming protogalaxy 
of intermediate (Steidel et al. 1999) or large mass 
(Granato et al. 2001). 

In cither model, the early mass distribution 
N(M,z) is directly related to N PS (M h ,z); the 
model (4B) yields the result represented by the 
thick solid line in fig. 2. 

4. QSs in interacting galaxies for z < 3 

After z w zq the galaxies are assembled into 
small groups of mass Mq 10 13 M Q , where the 
dominant member recurrently interacts with its 
companions, to the effect of refueling and rekin- 
dling an old BH; growing evidence (referenced in 
CV00) relates many QSs with interacting hosts. 
Small groups, with their high galaxies density and 
low velocity dispersion V still close to the galax- 
ian dispersion a, constitute preferred sites for such 
interactions to occur. 

Supply-limited accretion prevails here, since 
now the gas mass m in the host is depleted but 
no longer replenished by the interactions. Still, 
considerable fractions of the gas initially orbiting 
in the host at r ~ kpc are destabilized and partly 



made available for accretion; such fractions are 
easily evaluated (see CV00) in the form 

f d < Aj/j^GM'/aVb . (7) 

This includes the host structural parameter j/r, 
taken to be close to a; it also includes encounter 
orbital parameters: the impact parameter b, the 
relative velocity V, and the partner mass M'. 
Truly tidal interactions imply a postfactor r/b. 

With V a and b bounded by the group ra- 
dius, interactions and gas inflow both take times 
At w b/V <~ 10 _1 Gyr. The above equation easily 
yields fa in excess of a few %, of which about 1/3 
may reach the nucleus while the rest is likely to 
end up in circumnuclear starbursts, see Sanders & 
Mirabel (1996), CV00. This is enough to produce 
(see eq. 3) outputs L <; 10 46 erg s -1 in a host 
still gas rich with m <~ 10 10 M Q . Fractions up to 
fd ~ 1/2 are indicated by eq. (7), and in fact ob- 
tain in numerical simulations of grazing collisions, 
which arc statistically fewer. The full probabil- 
ity of accreting a fraction / is calculated by CV00 
to read P(f ) oc f~ 2 on the basis of eq. (7); this 
defines range and shape of the LF. 

In time, many small groups merge into richer 
ones where interactions are less frequent and ef- 
fective. Their volume density, grown rapidly later 
than zq, at low z goes down into a demise Nq{z) oc 
(1 + z); in addition, the decreasing interaction rate 
t 1 T 1 (z) « 0.5(1 + z) 3 / 2 Gyr _1 lowers the number 
of activated sources. The result for the QS pop- 
ulation is a moderate "density evolution" propor- 
tional to Ng(z)t~ (z). But a stronger "luminos- 
ity evolution" occurs since on average L oc fm(z) 
holds, and the residual host gas m(z) is depleted 
on a scale t ev as it is destabilized and used up in 
accretion episodes and by accompanying nuclear 
starbursts, with no replenishment. 

As L decreases and M increases, the sources 
on average go toward sub-Eddington luminosities, 
consistent with the data in Salucci et al. (1999), 
Wandel (1999). 

5. Relics at z ~ 

For early BHs growing at z > 3 the models (4A) 
or (4B) imply different scaling relations, that read 

Moc ( T 3 p- 1 / 2 (z)oc ( T 4 , (8 A) 
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or 

M cc a 5 , (8B) 

respectively. These stem from the simple hi- 
erarchical scaling of M/j oc V 2 R/G and R oc 
(Mh/p) 1 / 3 for the surrounding DM halos, as is ap- 
propriate for Mh < 10 13 M© when one galaxy per 
halo occurs. In such conditions, it is also fair to 
assume a oc V, and to fix the z-dependent fuzz ap- 
pearing only in eq. (8A) on relating the density to 

— 1/2 

Mh by means of p oc M h 1 , that holds for hierar- 
chically formed halos (see Haehnelt & Kauffmann 
2000). 

The relations (8A), (8B) turn out to be in tunc 
with the current debate concerning the MDO data. 
These have been recently recognized to follow a 
tight and steep correlation; the precise slope is 
given as slightly flatter than M oc cr 4 by Gcbhardt 
ct al. (2000), or close to M oc cr 5 by Ferrarese & 
Merritt (2000). In either case the scatter is found 
to be rather small, that is, factors 10 ±0 - 35 or less 
in M at given a. Note from fig. 1 that LFs flat as 
observed at high z obtain from model (4A) only 
upon convolution with scatter larger than 10 ±05 , 
as discussed by Haiman & Loeb (1998). 

Later than z « 3 fewer new BHs are still pro- 
duced, but the early ones still grow after M(z) = 
M(3) + J2i Ami, due to additional but dwindling 
accretion events Am^ caused by host interactions. 

The latter cause the mass distribution N(M, z) 
to evolve as given by 

d t N = / dfP(f) [N(M - fm) - N{M)\ . 

NtotTr J 

(9) 

Following up §4, Na{z)/N to t is the fraction of 
host galaxies residing in a group, relative to to- 
tal BH number including the dormant ones; for 
it we take values of order 10 _1 as discussed in 
CV00, based on the fraction 1/3 of bright galax- 
ies found by Ramclla ct al. (1999) to reside in 
groups with membership 3 or larger. The right 
hand side describes the net change at z of the 
distribution N(M, z), due to an accretion episode 
of Am = / m(z) which upgrades the initial mass 
M — Am to the current value M; the probability 
for this to occur is P(f), and the factor provided 
by the interaction frequency t~ 1 (z) converts it to 
a rate. 

The numerical solution plotted in fig. 2 shows 
how N(M, z) - starting from the condition at 



z = 3 computed in §3 - drifts and diffuses in time 
toward higher M, but only up to a cutoff; this is 
due to the large values of / being rare, and to m(z) 
being depleted from its initial value 10~ 1 Mf l - So 
the additions Am = / m(z) decrease on average, 
and large M are unlikely to grow much larger. 

We derive the corresponding M— a relation not- 
ing from cq. (7) that the gas masses available for 
accretion follow Am = jd m/3 oc m/a. On adopt- 
ing the scaling m(cr) oc cr 4 indicated by the Faber- 
Jackson relation or produced by stellar feedback, 
the result is Am oc cr 3 . If used in full, such masses 
dominate at low cr < 200 km s _1 , to yield maximal 
M values scaling as 

M oc cr 3 , (10) 

see fig. 3. The scatter is within an overall factor 5, 
the effective range given by the steep form of P(f). 
Truly tidal interactions (see the comment under 
eq. 7) yield a somewhat steeper scaling Am oc 
m((j) r I 'cr oc cr 4 and more scatter. 

But - given that enough gas is made available 
by the interactions as shown in §4 - the mass 
actually accreted may still be constrained by the 
QS feedback, depending on the degree of coupling 
of the source output to the surrounding baryons. 
Here we focus on the 90% radio quiet sources 
where the output is mostly radiative and roughly 
isotropic. At low z the unbinding constraint simi- 
lar to cq. (5) reads 

L At ^ nc 2 Am < e- 1 m{z) a 2 . (11) 

With e ~ 10~ 2 , the masses accreted in a host with 
a < 200 km s _1 are constrained to stay under 
those made available by interactions and given by 
eq. (7); using again m oc cr 4 , here the result is 
Ami oc cr 6 . These moderate additions drive con- 
siderably less evolution of N(M, z) as represented 
by the thin solid line in fig. 2; then M is always 
dominated by the initial values which follow eq. 
(8B), and will in fact converge to values scaling as 
M oc cr 5 . 

The full scaling laws given by eqs. (8B) with 
(11), or by (8A) with (10), are represented in fig. 
3. They differ mostly in the lower left section, 
where data are difficult to obtain from the kine- 
matics of stars or of a nuclear disk. 
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6. Conclusions and discussion 

We conclude that the relic BH masses are re- 
lated to the host galaxy dispersions by the steep 
correlation M oc a 5 if QS feedback onto the host's 
gas always provides the main constraint to the ac- 
cretion. In other words, this form is to hold at all 
a if the host potential well bounds the actual ac- 
cretion rates all the way from early BH formation 
to later additional growth. 

On the other hand, if such a control was never 
the limiting factor we expect the softer (and 
fuzzier) correlation M oc a 4 to hold at high a, and 
this to soften yet at lower a to the form Man 3 
set by the production rate of the available gas. 

These two extremes bracket the QS and BH 
story. The latter is centered on one kind of engine, 
the BH based on strong gravity; but it comprises 
different regimes of fueling, that is, of accretion 
and growth triggered by the environment under 
weak gravity. 

At early z ?J 3 much fresh gas is made available 
for full, self-limiting accretion by major merging 
events which build up the halos and the embedded 
spheroids. BHs grow rapidly, generating the up- 
per section of the correlation M oc a 4 -f- a 5 . Mean- 
while, the associated QSs flare up at Eddington 
luminosities; their LF grows in range and height, 
tracking the progressive halo build up. 

Later than z w zq ~ 2.5 the dearth of gas curbs 
the accretion and bends down the QS evolution, as 
host interactions with group companions are still 
able to trigger accretion but no longer to import 
fresh gas. Though supply-limited, the gas masses 
made available for accretion are still sizeable. If 
these are used in full, the average luminosities are 
still high and the BHs still grow considerably; re- 
latedly, the correlation is softened to M oc a 3 in 
its lower section. 

At z lower yet, the host gas reservoirs approach 
exhaustion, and weaker AGN activity may be sus- 
tained by lesser gas productions or smaller sup- 
plies provided, e.g., by internal instabilities (see 
Heller & Shlosman 1994, Merritt 1999) or by satel- 
lite galaxies cannibalized by the hosts (CVOO). 
These latest and smaller accretion events can fuel 
many weaker AGNs more likely to be pinpointed 
in X-rays, implying widespread accretion power 
but moderate growth of individual BHs, as also 
argued by Salucci et al. (1999). Deep surveys 



in hard X-rays are now uncovering optically hid- 
den, and even X-ray obscured accretion power, as 
they resolve and count the AGNs that dominate 
the hard XRB (see Giacconi et al. 2001, Hasingcr 
ct al. 2001); the indications are currently favor- 
ing AGNs of intermediate powers and redshifts as 
main XRB contributors (see Comastri et al. 2001). 
These findings jointly with the scanty evidence 
of, and the indirect bounds to a QS 2 population 
(see Maiolino et al. 2001) suggest to us that the 
optical-IR sampling is fair to within a factor 2 as 
for the outputs of early and luminous QSs, and is 
reliable for linking these with the largest BHs. 

Thus our conclusions may be rephrased in 
terms of two simple patterns standing out of this 
complex picture. 

First, the largest MDOs in local inactive galax- 
ies are directly linked with the high-z luminous 
QSs. In particular, we have shown that flat opti- 
cal LFs given by the feedback-constrained model 
(4B) arc linked to the steep and tight correlation 
given at large masses by M oc a 5 (see cq. 8B). 
Current evidence concerning high brightness of the 
hosts (see §3), flat shape of the high-z LFs (§3), 
and narrow scatter in the M — a correlation (§5) 
concur to favor this specific link. But telling data 
will be provided by the precise slope of the upper 
MDO correlation observed locally in the optical 
band, jointly with the rate of QS decline toward 
high z derived from large area surveys out to the 
near-IR like the Sloan Digital Sky Survey. 

Second and concerning the lower region of the 
M-a plane, it is conceivable (within the present 
uncertainties about e) that the fuel throttle to the 
BH engine is always controlled by the host wells as 
is expressed by eq. (11); if so, the main accretion 
modes may be effectively united by the feedback 
action. In §5 we have shown that this provides the 
lower bound to M expressed again by M oc a 5 . 
On the other hand, an upper bound M oc a 3 is 
set by the condition of maximal fuel availability, 
eq. (10); this includes the production of optically 
obscured accretion power. Note that the highly 
obscured, wind blowing BHs suggested by Fabian 
(1999) would be located in between. The actual 
feedback action will be observationally probed by 
the slope of the M — a correlation in its lower sec- 
tion; here the growing reverberation map database 
(see Ferrarese et al. 2001) in current AGNs will 
be telling. 
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Fig. 1. — High-z LFs from the models discussed in 
§3: eq. (4A), dotted; eq. (4B), solid. Here Lb is 
the luminosity in the blue band (bolometric cor- 
rection k = 10), and L* = 10 45 erg s _1 ; canonical 
ACDM cosmology/cosmogony. Data: filled circles 
from Kcnnefick et al. 1995, Schmidt et al. 1995; 
open symbols from Fan et al. 2001a. 
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Fig. 2. — The evolution of the mass distribution 
N(M, z) of the relic BHs. Down to z w 3 this 
is given by Nps(Mh, z) with the transformation 
(4B), see the thick solid line. Afterwards, the 
distribution is driven by host interactions as de- 
scribed by eq. (9). The dashed line corresponds 
to unconstrained accretion of 1/3 of the gas sent 
toward the nucleus; the thin solid line corresponds 
to accretion constrained by the QS feedback fol- 
lowing eq. (11). The related local mass density 
amounts to pbh ~ 4 and 210 14 M Q Gpc~ 3 , re- 
spectively. Cosmology/cosmogony as in fig. 1. 
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Fig. 3. — The models for the local M — a relation 
discussed in §5. Dashed: unconstrained accretion 
at levels / = 5%, the average value given by the 
distribution P(f). Solid: feedback constrained ac- 
cretion with e = 5 10~ 3 . Data from Gebhardt et 
al. 2000, and Ferrarese & Merritt 2000, with the 
open symbols marking results from reverberation 
mapping. The model normalizations at high M 
follow the fits by the same authors. 



